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ABSTRACT 

Magnetohydrodynamic accretion disk simulations suggest that much of the energy liberated by the 
magnetorotational instability (MRI) can be channeled into large-scale toroidal magnetic fields through 
dynamo action. Under certain conditions, this field can dominate over gas and radiation pressure in 
providing vertical support against gravity, even close to the midplane. Using a simple model for the 
creation of this field, its buoyant rise, and its coupling to the gas, we show how disks could be driven 
into this magnetically dominated state and deduce the resulting vertical pressure and density profiles. 
Applying an established criterion for MRI to operate in the presence of a toroidal field, we show that 
magnetically supported disks can have two distinct MRI-active regions, separated by a “dead zone” 
where local MRI is suppressed, but where magnetic energy continues to flow upward from the dynamo 
region below. We suggest that the relative strengths of the MRI zones, and the local poloidal flux, 
determine the spectral states of X-ray binaries. Specifically, “intermediate” and “hard” accretion 
states occur when MRI is triggered in the hot, upper zone of the corona, while disks in “soft” states 
do not develop the upper MRI zone. We discuss the conditions under which various transitions should 
take place and speculate on the relationship of dynamo activity to the various types of quasi-periodic 
oscillations that sometimes appear in the hard spectral components. The model also explains why 
luminous accretion disks in the “soft” state show no signs of the thermal/viscous instability predicted 
by standard a—models. 

Subject headings: accretion, accretion disks — binaries: close — black hole physics — magnetic fields 
— X-rays: binaries 


1. INTRODUCTION 

Despite many successes, standard accretion disk mod¬ 
els fail to account for key features of X-ray binary 
(XRB) phenomenology. Both analytic a—viscosity mod¬ 
els (jShakura fc Sunvaevlll97.‘f l and numerical simulations 
capturing angular mom entum transport by ma gnetoro- 
tational instabili ty le.g.lB albus fc Ha.wlevH l 9981 ) predict 
strong thermal (|Jiang et all l2014bh and viscous insta¬ 
bility when radiation pressure dominates other forms of 
pressure in the disk, yet no evidence of such instability 
is obviously seen during XRB outbursts. Standard mod¬ 
els also fail to provide a satisfactory explanation for why 
accretion flows undergo spectral state transitions, where 
the spectrum switches from soft and thermal to hard 
and nonthermal, and vice-versa. Models for the transi¬ 
tion typically rely on the existence of t wo possib l e disk 
structures at th e same accretion rate dChen et all 119951 : 
[Esin et MHT997 ) — one cool, dense and thermal and the 
other hot, tenuous and out of equilibrium — but provide 
no compelling mechanism to explain why one is chosen 
over the other when conditions are such that either could 
exist. This problem is most acute at soft-to-hard tran¬ 
sitions, which occur at relatively low luminosities where 
thin disks are expected to be relatively stable. One sug¬ 
gestion to explain this transition posits evaporation of 
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the disk by a hot corona dMev er fe Mever-Hof meisted 
11994 iMever et all 1200(1 : iRozariska. fc Gzernvl 12000( 1 . but, 
begs the question of what keeps the corona heated at 
such a low density that electrons and ions are thermally 
decoupled, which appears to be necessary in order to pre¬ 
vent the electrons from radiating away most of the heat. 
Another difficulty is that the states are not mutually ex¬ 
clusive, but often seem to co-exist in relatively long-lived, 
hybrid “intermediate” states. 

Non-standard disk models that include magnetically 
dominated regions provide a possible solution to these 
problems. At a minimum the buoyant rise of magnetic 
held generated within the disk — which in weakly mag- 
netized disks is typically a cyclic, dynamo - like process 
((Brandenburg et all I 1 9951 iStone et all 1 1 9961 iDavis et~ahl 

l2010HQ’Neill et al.ll2011l : lSimon et al .1120121) — can trans¬ 

port substantial amounts of energy vertically, and could 
heat low-density gas to create a corona. Turbulent ver¬ 
tical energy transport, moreo ver, has a stabilizing influ¬ 
ence on thermal instability dZlru fc Naravanl 12014 ). It 
is at best unclear, however, whether a disk model that 
sandwiches a gas or radiation-pressure dominated mid- 
plane with a magnetically dominated corona can undergo 
a state transition. Fully converting a thin, optically thick 
accretion disk to a tenous hot state requires such a dra¬ 
matic decrease in density that it is hard to see how it 
could be accomplished without some catastrophic event 
to destroy the disk, and such events appear incompati¬ 
ble with observations t hat suggest rather seamless tran¬ 
sition s among states (|Dunn et al.1 l2010fc iKalemci et al.1 

HhJ). 

In this paper, we suggest that the full catalog of 
XRB accretion states can be understood in terms of a 
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competition between the dynamo generation of toroidal 
field and its buoyant rise, provided that the disk mid¬ 
plane is magnetically dominated. Magnetically dom¬ 
inated disk solutions can be constructed in a man¬ 
ner that parall els the s tanda rd Shakura-Sunyaev (1973) 
models (IPariev et al.1 120031b though it is not obvi- 
ous analy tically jwhat sets the level of magnetization. 
iBeeelman fe Pringle! (l2007f) presented models for max¬ 
imally magnetized, magnetically supported disks, in the 
limit where radiation pressure dominates over gas pres¬ 
sure. The maximal level of magnetization was cho¬ 
sen to satisfy the marginal criterion for axisymmetric 
MRI to operate in a disk with a suprathermal toroidal 
magnetic field (j>b > Pg as) and a weak poloidal field, 
v a < c g v K, where v\ is the Alfven speed associated 
with the toroidal field, c g is the sound speed associated 
with the gas pressure on ly, and vk is the Keplerian speed 
(jPessah fc Psaltisl 120 05). Magnetically dominated disks 
are known to have various phenomenologically interest¬ 
ing properties. Th ey are stab le agains t ther mal and vis¬ 
cous instability (IBeeelman fc Pringle! 120071 ). apparently 
persistent if they can b e established in the first place 
( Johanseii^JjCvin 2008 j) , and resistant to fragm entation 
( Begelman fc Pringlell2007t iGaburov et al 1120121 b 

Numerical simulations show that a magnetically dom¬ 
inated disk forms w henever th e p oloidal flux exceeds a 
threshold value (|Bai fc Stonell2013lf , given by the require¬ 
ment that the pressure in the poloidal field reaches about 
a percent of the sum of the gas and radiation pressures. 
When this condition is met, the toroidal field becomes 
the dominant form of pressure throughout the disk, sup¬ 
porting the gas against the vertical component of gravity. 
Here we study the implications of the resultant vertical 
structure for XRB states. The key result is that magnet¬ 
ically dominated XRB disks may possess a single MRI- 
active layer, centered on the equator, o r two active la yers 
separa ted by a dead zone in which the iPessah fc Psaltisl 
(120051) criterion for MRI activity is violated. When a 
dead zone exists, substantial amounts of magnetic en¬ 
ergy can be injected into the gas at the very low den¬ 
sities where electrons and ions are thermally decoupled, 
producing a spectrally hard flux component that can co¬ 
exist with a cool equatorial region over a range of ac¬ 
cretion rates. In its simplest form, the state of the disk 
is a local function of the accretion rate and a, which 
is itself determine d by t he streng th of the net poloidal 
flux. Following [Begelman fc Armitagel (|2014f) . we appeal 
to global changes in the net flux to explain the hysteresis 
seen in transitions between the states. 

In subsequent sections, we present a simple analytical 
model for the vertical structure of accretion disk+corona 
systems, and apply this stratified model to XRBs. In 
Section 2 we show that the competition between toroidal 
field generation and escape is finely balanced, with small 
changes in dynamo and buoyancy parameters leading 
to large changes in the ratio of magnetic pressure to 
gas+radiation pressure near the equator. We estimate 
the gradients of magnetic pressure and density, both 
when MRI is active and when it is suppressed. In 
Section 3 we use this model to extract the basic phe¬ 
nomenology of magnetically supported disks, includ¬ 
ing the likely suppressi on of MRI in regions where the 
IPessah fc Psaltisl (|2005H criterion is violated. We ap¬ 


ply these results to XRBs in Section 4, where we ar¬ 
gue that the relative amounts of energy released in the 
lower and upper MRI zones determine whether the flow 
is in a soft, int e rmedi ate or hard state. (We note that 
iRozariska et al.l (1201511 arrive at a similar vertical struc¬ 
ture, that includes a dead zone, from spectral fitting.) By 
linking the vertical structure to the outcome of dynamo 
activity in magnetically supported disks, our model pro¬ 
vides an explanation for the rich phenomenology of XRB 
state transitions. We discuss the results and conclude in 
Section 5. 

2. TOROIDAL FIELD GENERATION AND ESCAPE 

Consider a steady-state disk atmosphere at cylindri¬ 
cal radius R in the gravitational field of a central mass 
M, supported vertically by some combination of ther- 
mal+radiation pressure, p(z), and magnetic pressure due 
to a toroidal field, ps(z) = B^/8n. Denoting the vertical 

gravity by g z « GMz/R 3 = fl 2 z, we have the equation 
of hydrostatic equilibrium, 

■^(p + Pb) = -ptfz, ( 1 ) 

where p is the density. A second equation is provided by 
considering the vertical flux of magnetic energy (Poynt- 
ing flux), which we treat phenomenologically according 
to a set of physically motivated assumptions about its 
creation and escape. Specifically, we suppose that MRI 
liberates accretion energy locally at a rate (per unit vol¬ 
ume) ol£1(p+pb), of which a fraction as/ct is channeled 
into the toroidal field. This toroidal field escapes, via 
the Parker instability, with a mean vertical speed tb(-z), 
and thus carries a mean Poynting flux Fb{z) = 2pbVb- 
Because the configuration is expected to be far from 
marginal stability, we expect the rise of field lines to re¬ 
flect the near free-fall of gas sliding along buckled field 
lines, and thus take 

vb{z) = 7i£Iz, ( 2 ) 

where f Iz is the local free-fall speed at 2 and p < 1 ac¬ 
counts for geometric effects, such as field-line tangling, 
finite time scale of reconnection, etc., that might impede 
the overall energy transport. We assume (at least ini¬ 
tially) that r) is independent of z but it need not be a 
universal constant. 

Rather than lifting the gas uniformly, the field repeat¬ 
edly buckles (probably on a transverse scale larger than 
the local height: Tout & Pringle 1992), allowing mat¬ 
ter to slide downward as it rises. Large-scale toroidal 
loops of field then reassemble themselves slightly higher 
in the disk atmosphere, through reconnection, and the 
process repeats. This process both slows down the net 
rise of the field, compared to the instability time scale, 
and causes the magnetic field to lose energy to the gas as 
it rises, at a rate (per unit volume) —VBdpB/dz. Some of 
this energy goes initially into potential energy of the gas, 
then kinetic energy as the gas drops back, and ultimately 
internal energy as the kinetic energy is dissipated. Addi¬ 
tionally, magnetic energy can be converted directed into 
kinetic energy and heat at the reconnection sites. We 
gloss over these details and assume that all of the en¬ 
ergy lost by the Poynting flux goes into heating the gas 
locally. An analogous heating prescription is frequently 
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used in treatments of cosmic-ray transport in which cos¬ 
mic rays stream through the background plasma but are 
isotropized in a frame co-moving with an ensem ble of 
Alfven waves that scatter the cosmic rays ljSkillindll975j l. 
In this case, the cosmic-ray streaming speed (basically a 
weighted average wave speed) plays the role of vb- An 
equivalent description was used to describe “effervescent 
heating” of the intracluster medium of a cluster of galax¬ 
ies by buoyant bub bles of hot p l asma i njected by a cen tral 
radio galaxy fBcgelman 2001; Ruszkow ski fe Beeelmanl 
12001 . 

If there are no additional energy losses, we may write 
down the following equation for the Poynting flux, 


cLFb 

dz 


dp B , „ dv B 

-E 2 nn- 


d,7)V! 


(3) 


equation of state or a radiative diffusion equation. How¬ 
ever, some generic features of the pressure profile are in¬ 
sensitive to the detailed behavior of the density and can 
be gleaned from the asymptotic behavior of solutions at 
small and large heights above the equator. For example, 
by assuming that p = 1 + 0(y) for y <C 1, we find that 
the normalized pressures close to the equator are given 

by 


Pb 


1 - 


Po 

2(1 + A)) y ’ 


p 


1 - 


(1 + 2/3o) 
2(1 + A,) y 


( 10 ) 


This implies that Pb — p ~ y/ 2 increases with height, 
at least close to the equator. In fact, we can prove that 
Pb/p > 1 at all heights, by combining equations ([ 6 ]) and 
© to yield 


which simplifies to 


-JZ^PbVb) = a B ttv B {p + Pb)- (4) 

Defining the central values of pressure and density, 
Po = p(0), Peso = Pb( 0) = A/Vo, P° = p( 0 )j a § as Pres¬ 
sure scale height H = f ^~ 1 (po/Po) 1 ^ 2 , and a dimension¬ 
less vertical variable 


V = 


2H 2 1 


(5) 


and using eq. © for ub, we obtain the dimensionless 
equations of hydrostatic equilibrium and magnetic en¬ 
ergy, 

-t-(Pb + A )P) = ~Po P (6) 

ay 

-r-(PB2/) = ^-{pb+PoP), (7) 

ay 277 


where a tilde denotes a quantity normalized to its equa¬ 
torial value, so that pb(0) = p( 0 ) = p( 0 ) = 1 . 

Since equation © is singular at y = 0 , the parameters 
must satisfy a regularity condition to ensure that pg( 0 ) 
is finite (where a prime indicates differentiation with re¬ 
spect to y), 

A) = — - I- (8) 

a b 

We can therefore write eq. © as 

VPb = 1 ip - Pb)- (9) 


Equation j5|) reveals an important generic feature of the 
model: the ratio of gas+radiation pressure to magnetic 
pressure on the equator, /Jo, depends solely on the ratio 
of the escape speed parameter p to the flux production 
parameter as- When this ratio is 7§> 1, flux can read¬ 
ily escape as it is produced and gas+radiation pressure 
dominates the disk. When this ratio equals 1, the pres¬ 
sures are equal, but only a small decrease in the ratio, to 
1 / 2 , drives the ratio of gas+radiation pressure to mag¬ 
netic pressure formally to zero. Thus, an order unity 
change in a simple parameter ratio makes an enormous 
difference to the physical nature of the disk. 

To determine the pressures as a function of height re¬ 
quires a third equation involving the density, such as an 


(1 + Po) 


dhi(p B /p) 

dlny 





( 11 ) 


Since the right-hand side of eq. CD is positive definite 
for all pb < P, we see that the ratio always tends to 
values > 1 . 

If we assume that pb ^ P at large y, we obtain the 
asymptotic scaling Pb oc y~@ 0 /O-+P 0 ) ^ w j 1 i c i 1 we can com¬ 
bine with eq. m to obtain a single expression that in¬ 
terpolates between the two limits, 


/ y \ ~ /W(l+/3o) 

PB « (l + |) . (12) 

For a polytropic equation of state, p = p 7 , eq. (fl2l) (which 
does not depend on 7 ) gives an extremely good fit to the 
normalization as well as the slope of pb at large y (Fig. 1), 
implying that the magnetic pressure is insensitive to the 
gas equation of state. Not surprisingly, the gas pressure, 
and thus the density, does depend on the equation of 
state (Fig. 2). 

For polytropic equations of state, magnetic pressure 
dominates above a few scale heights (Fig. 2), even when 
gas+radiation pressure is dominant on the equator. In 
fact, for these models the “temperature” p/p decreases 
with height, whereas the effective virial temperature in¬ 
creases oc z 2 . Thus, these coronal models are “cold” and 
fully supported by the gradient of magnetic pressure. 

A remarkable feature of these models is that relative 
importance of coronal vs. disk heating depends entirely 
on the ratio of gas+radiation pressure to magnetic pres¬ 
sure on the equator, Po- The vertically integrated heat¬ 
ing rate from the equator up to z H converges to a 
constant for a gas+radiation dominated disk, Po > 1 , 
whereas it continues to increase oc PbZ oc ,A 1- / 3 °)/( 1 +/ 3 °) 
for Po < 1. Thus, in this model the sensitive switch, 
eq. ©, that determines Po , also determines whether most 
of the dynamo energy is dissipated below the disk photo¬ 
sphere, or travels far above the disk plane into a region 
where optical depths are low and cooling may be ineffi¬ 
cient. 

Frequent reversals of the toroidal field, which are ob- 
serve d in accretion disk simulations with a net magnetic 
flux (IDavis et all I 2010 L IQ ’Neill et all 120111 : ISimon et all 
l 2012 f) can steepen the magnetic pressure gradient at large 
z, due to increased magnetic reconnection across the ris¬ 
ing current sheets that separate the reversed regions. In 
an Appendix, we show how to incorporate such rever- 
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Figure 1. Normalized mag netic pressure pe from solutions of equations (( 6 } and J9} for an equation of state p = fp, compared with the 
analytic approximation G2J- Ordering of solid curves (top to bottom, along right-hand side) corresponds to 7 = 1,4/3, 5/3. Dashed curve 
is analytic approximation. Left: Gas+radiation pressure-dominated case, j3 0 = 10. Right: Magnetic pressure-dominated case, (3q = 0.1. 



H H 

Figure 2. Gas+radiation pressure (solid lines) relative to normalized magnetic pressure (analytic approximation G3: dashed lines) as 
a function of height, for an equation of state p = p 7 . Left: Case where gas+radiation pressure dominates on the equator, /3q = 10. 
Right: Magnetic pressure-dominated case, (3o = 0.1. Ordering of solid curves (top to bottom) corresponds to 7 = 1,4/3, 5/3. Even when 
gas+radiation pressure dominates on the equator, magnetic pressure becomes dominant within a few scale heights. 


sals into our model. However, the dynamos producing 
magnetically supporte d (0g < 1) disks ap pear to reverse 
infrequently, if at all, dBai fe Stonell2013l Salvesen et al., 
in preparation), allowing us to apply models with a uni¬ 
directional field. 

3. PHENOMENOLOGY OF MAGNETICALLY SUPPORTED 

DISKS 

In this section we extract what we believe are the 
salient lessons about the generic behavior of magneti¬ 
cally supported disks and coronae, which will form a 
basis for our applications to X-ray binary state transi¬ 
tions. Since we are concerned here with magnetically 
supported disks that exhibit infrequent toroidal field re¬ 
versals dBai fc StonellMl . we will ignore the effects of 
reconnection losses across large-scale current sheets. 

3.1. When are disks magnetically dominated? 

We have shown that the degree of magnetic domination 
can be described by the difference between two parame¬ 
ters, 2?;—ttB, where ap describes the rate of field creation 
and r/ describes its rate of loss through buoyancy. Ob¬ 
viously these parameters cannot be physically indepen¬ 
dent, since if ap approached 2r] the magnetic pressure 
would become infinitely strong compared to other forms 
of pressure. Before this limit could be reached, addi¬ 


tional field would be expelled, or the dynamo quenched, 
leaving 2?? at least slightly larger than ap. 

Simulations suggest that the difference between these 
parameters is governed by the poloidal 0 —parameter 0 p , 
i.e., the ratio of the gas+radiation pressure to the pres¬ 
sure in the poloidal components of the magnetic field 
(Bai & Stone 2013, Salvesen et ah, in preparation). 
This parameter is also thought to govern the effective 
value of the viscous a—parameter, with a reaching val¬ 
ues ~ 0(1) for j3 p 10 2 , while settling at much lower 
values ~ 0.03 for 0 p 10 4 . Similarly, the toroidal 
magnetic pressure associated with the dynamo is very 
small for large values of 0 P , begins to approach the 
gas+radiation pressure near the equator for 0 p 10 3 , 
and becomes dominant throughout the disk (by about a 
factor ~ 10) for 0 p ~ 10 2 . These poloidal field strengths 
fall squarely in the regime where local numerical simu¬ 
lations pr edict that magneto-centrifugal winds would be 
l aunc hed (|Blandford fc Pavnel 119821 : iSuzuki & Inutsukal 
I2009T ). though the local approximation is ill-suited to 
determini ng the resu lt ant mass and angular momentum 
loss rates ijFromang et al.ll2014f) . Wind-driven flux trans¬ 
port would assist in the maintenance of the relatively 
strong poloidal fields needed for our model, though the 
details of such transport would need a global treatment 
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that we do not attempt here. 

For the remainder of this paper we will assume that 
the strength of the dynamo, hence the dominance of the 
toroidal magnetic field, is set by some external parame¬ 
ter, presumably /3 p . We will characterize the dominance 
of the toroidal pressure in terms of an equatorial value 
of the toroidal /?—parameter, ft 0 . 

3.2. Approximate vertical structure 

In regions where MRI is active and the disk is mag¬ 
netically supported, we will assume that the magnetic 
pressure above 2 ~ H varies as a weak power of z, 
Pb oc z~ 2 ^ a , and therefore can be treated roughly as 
a constant for ft'Cf Maintenance of approximate hy¬ 
drostatic equilibrium then demands that p oc z~ 2 . 

Where MRI is suppressed, the gradient of magnetic 
pressure is governed by losses in the vertical Poynting 
flux, according to eq. d4]) with «b = 0. We then have 
Pb oc z~ 2 and p oc 2 ~ 4 . 

In such “dead zones,” the Alfven speed tracks the local 
free-fall speed, v\ = 2 pb/p = Q 2 z 2 . In MRI-active re¬ 
gions, however, the weak pressure gradient implies that 
the local Alfven speed exceeds the free-fall speed by a 
factor ~ Pq 1 . Taken at face value, these scalings would 
imply a sharp density inversion at the transition from an 
MRI-active zone to a dead zone, but we expect such an 
inversion to be unstable, and to be smoothed, e.g., by 
a local increase in the buoyant rise speed or a gradual 
rather than sudden drop in ae- For purposes of outlin¬ 
ing our state transition model, however, we will assume 
the simple piecewise scalings outlined above. 

3.3. Suppression of MRI 

In order for MRI to operate in the presence of both 
poloidal and toroidal magnetic field, two conditions must 
be satisfied. First, the poloidal magnetic pressure must 
be smaller than the pressure responsible for setting the 
scale height, which in this case means that the poloidal 
field must be smaller than the toroidal field. We assume 
that this is always satisfied near the disk midplane, since 
the presence of a moderate poloidal field appears to be 
capable of catalyzing the creation of a much stronger 
toroidal field through dynamo action. At large heights, 
however, this condition might not be met if the toroidal 
magnetic pressure drops significantly across a dead zone, 
since poloidal field strength should be roughly indepen¬ 
dent of z. 

The second condition provides an upper limit to the 
toroidal field. Pessah & Psaltis (2005) showed that MRI 
is suppressed, in the presence of a strong toroidal field 
and weak poloidal field, when 

v\ > c g v K (13) 

where c g is the thermal (gas) sound speed and vk is 
the Keplerian speed. Begelman & Pringle (2007) inter¬ 
preted this criterion in terms of accretion disk structure, 
and argued that the sound speed should correspond to 
that of the gas only and should not include contributions 
from radiation pressure because of the highly diffusive 
nature of radiative transport across small scales. They 
also noted that this criterion is consistent with earlier 
resul ts o n no n- axisymmetric MRI for a purely toroidal 
field (iTerauem fe Papaloizoulll996l f. We are not aware of 


rapid accretion through 
two temperature zone 



Figure 3. Illustration of the vertical structure of a magnetically 
dominated disk in the regime which we identify with the interme¬ 
diate state of XRBs. A strongly magnetized MRI-active mid-plane 
is separated from an MRI-active two temperature corona by a dead 
zone where the toroidal magnetic field is too strong to admit in¬ 
stability. 

any numerical simulations that have directly tested the 
Pessah fe Psaltis (2005) criterion, though the results of 
IGaburov et all (120121) show that magnetically dominated 
disks are both relatively stable an d can support MRI in 
the supra-thermal disk core. The IGaburov et all (120121 1 
simulations, however, did not include the extra physics 
(such as the possibility of separate electron and ion tem¬ 
peratures) that we argue is essential for a model of XRB 
state transitions. 

4. STATE TRANSITIONS IN X-RAY BINARIES 

4.1. Schematic approach to state transitions 

We consider a magnetically supported disk-corona sys¬ 
tem that is MRI-active close to the equator, and follow 
its behavior with increasing height. Since the v\ ~ z 2 , 
eq. m will be satisfied above some height zi, shutting 
off MRI, unless c g oc T 1 / 2 increases very rapidly with 
height, where T is the temperature. In the optically 
thick, thermalized equatorial region of an accretion disk, 
T is a relatively weak function of height, and therefore 
the system enters a dead (MRI-inactive) zone at 2 > Z\. 

In the dead zone, gas continues to be heated by the 
rising toroidal flux from below, but at a rate that de¬ 
creases with height ~ pb oc 2 -2 . v\ continues to increase 
oc 2 2 (albeit with a slightly different normalization), but 
the most dramatic change is the steep decline in density, 
p oc 2 -4 . For the case of X-ray binaries (XRBs), we will 
argue that the decline in the cooling rate, due to the 
decrease in density, far outweighs the weakening heating 
rate and leads to the thermal decoupling of electrons and 
ions and a runaway increase of ion temperature, 7] T e . 
Remarkably, we find that c g oc T ; ' can increase faster 
than 2 2 , leading to the re-establishment of MRI above 
a height 22, provided that the poloidal magnetic field is 
not too strong. We note that the energy source for this 
second MRI zone is not the Poynting flux generated near 
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the equator, but rather the local differential rotation of 
the flow within the zone, and the gravitational binding 
energy liberated by accretion of this gas. However, the 
Poynting flux from below is crucial for catalyzing the re¬ 
establishment of MRI by heating the gas sufficiently to 
drive it back into the MRI-unstable regime. Once estab¬ 
lished, this upper MRI zone extends all the way to the 
top of the corona at z ~ R. 

Because pb is roughly constant with height within each 
MRI-active zone, the total dissipation within that zone is 
~ cxQpb (^max - 2 m in), where z m in and z max are the lower 
and upper boundaries. For simplicity in the approximate 
expressions derived below, we will assume z max z m 
Then, given that an upper MRI zone exists, Z 2 < R, 
we can estimate the ratios of accretion rates and powers 
associated with the two MRI-active zones: 

M2 _ L-2 ~ PB 2 R ^ Z\R 

Mi Li pbiZi z\' 

where the final relation reflects our extremely crude ap¬ 
proximation that pbi!pb\ ~ ( 21 / 22 ) 2 . In this picture 
the two MRI-active zones share the same poloidal flux, 
but are otherwise independent. While we remain in a 
regime where the dominant toroidal field controls MRI- 
activity, this is locally a reasonable assumption. Glob¬ 
ally, however, neighboring regions of the disk will be cou¬ 
pled by the evolution of the poloidal field, which will be 
qualitatively different depending on whether the verti¬ 
cally _extended upper MRI-active layer is present or not 
dLovelace et all 120091 ). We expect inward transport of 
the net poloidal flux to be greatly facilitated when there 
is a fast-accreting upper layer, though the details will de¬ 
pend substantially on the strength and nature of MHD 
outflows which this geometry also favors. 

Schematically, we associate the low/hard, interme¬ 
diate/hard, intermediate/soft and high/soft accretion 
states with the cases M 2 /Mi 1, ;> 1, 1, and 0, 

respectively. In the following section we estimate these 
criteria quantitatively for conditions appropriate to X- 
ray binaries, and flesh out some of the associated obser¬ 
vational characteristics. 


where 

% ~ Oma-^rhl^m-^x 1 / 36 ( 15 ) 

R 

(iBegelman fc Pringlel 12007 1. 

This result is based on the assumption of L TE ( the 
valid ity of which is discussed in IBegelman fe Pringle! 
(|2007ll 1. including thermal equilibrium between electrons 
and ions, T « T e . We can check the latter assumption 
by using a simple approximate expression for the volume 
energy transfer rate from ions to electrons, minimally 
given by Coulomb collisions: 




36 



5/2 


p 2 c 3 n 



/ Tj \ 1/2 T; — T e 
\TJ Te 


( 16 ) 


(jGuilbert fc Stemievl Il985f ). where we have taken the 
Gaunt factor In A in their original expression to be ss 25. 

Assuming that much of the direct MRI heating goes 
into ions and equating aflpB to Hi e , we obtain 


A = ^ 1/2 h-h „ 0.03a 29 / 18 m- 2/9 m- 7/1 V 47 / 36 , 

(17) 

implying that the species are well-coupled at 2 ^ Z\. 

At z ^ zi, the magnetic pressure goes from being 
roughly constant to declining oc z~ 2 , while the density 
tends toward a steep decline oc z~ 4 . Heating by MRI 
is gone, but the dissipation of Poynting flux heats the 
gas at a rate — vb( dps/dz) oc pb- Denoting conditions at 
Z\ by a subscript 1, assuming that the changes in den¬ 
sity and pressure slopes are instantaneous, and ignor¬ 
ing changes in normalization of the pressure and density 
power-laws, so that pb ~ pbi{z/zi)~ 2 , etc., we obtain 
A(s) ~ (Cg/Cgi)(z/zi) 6 Ai, implying that thermal cou¬ 
pling breaks down very quickly, i.e., at 2 <; 2z\ at this 
crude level of approximation. 

Above the decoupling point but still in the dead zone, 
we can assume Tj T e and write c g /c ~ (AT)/ pc 2 ) 1 / 2 , 
where p ss 0.6m p is the mass per particle. Writing 0 e = 
kT e /m e c 2 , we can express the thermal balance for ions 
in the form 


4.2. Soft state 

The structure of thin, magnetically supported disks 
with radiation pressur e exceeding gas pressure is dis¬ 
cussed at length in IBegelman fc Pringlel (j2007H . under 
the assumption that the disk marginal l y satis fies the in¬ 
stability criterion of IPessah fc Psaltisl (I2005f ). Readers 
are referred to that paper for insight into results pre¬ 
sented in this section. As we saw earlier, the vertical 
structure of a magnetically dominated disk guarantees 
that most of the accretion will occur close to marginal 
(maximally magnetized) conditions for MRI, although 
the optically thick disk interior may hide a large column 
density of slowly accreting matter. We denote the di¬ 
mensionless accretion rate through the equatorial layer 
as mi = Mi/MEdd, where M^dd = -bEdd/c 2 = 47 tGM/kc 
is the Eddington accretion rate, without the often-used 
efficiency factor. Here M = IOtoiMq is the mass of the 
central object and k is the opacity. We scale the radius 
to R g = GM/c 2 , x = R/Rg. The upper boundary of 
thi s accretion layer, where MRI shuts off according to 
the IPessah fc Psaltisl (I2005H criterion, is then given by zi 


^ - 0.01 6 3 J 4 (j-J* a 31 / 36 m- 2/9 m- 5/3 V 31 / 72 . (18) 

Now consider the criterion for the onset of MRI, v\ < 
CgVK- If we make the reasonable assumption that T e is 
increasing with z, then the right-hand side of this rela¬ 
tion increases at least as rapidly as z 3 . But the quantity 
on the left-hand side is constrained, by the equation of 
hydrostatic equilibrium, to rise oc z 2 . Therefore, there 
is a chance that MRI will be re-established, above the 
dead zone, under conditions very much hotter and more 
tenuous than in the equatorial accretion zone. 

Denoting Z 2 as the height where MRI is reestablished, 
we obtain 

/ a \ -3/4 

Z .l „ 0.4 ( AA ) a- 31 /36 rh 2/9 m 5/36 ;c 85/72_ (ig) 

The quantity d e \ depends on the temperature at z 1 , 
where electrons and ions are still coupled. A number 
of factors may affect the value of 9 e at 22 , but to ob¬ 
tain a very rough estimate, we assume that the electrons 
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close to Z 2 are cooling by unsaturated Comptonization 
so that the Compton y —parameter <, 1, depending on 
the amount of cooling required. If we further assume 
that the electrons are subrelativistic and the scattering 
optical depth r > 1, then 



y ~ 4t 2 6» 6 ~ 4r 2 6» e , 

(20) 

where 

ti ~ 300a -8 / 9 TO 7 / 9 my 9 x -5 / 9 

(21) 

(Begelman & Pringle 20071). We obtain 


9 e 2 r 

-0.03y 2 / 11 a- 73 / 99 77 l f /99 TO- 1/99 x 145 / 198 

(22) 

and 



Zl 

^5y- 3 / 22 a- 83 / 198 TOf /99 TOl /198 x- 25 / 396 . 

(23) 


Note the extreme insensitivity of this expression to all 
parameters except a and to — in particular, this result 
hardly depends on the value of y , which is uncertain. We 
also have 

^ - 0.3a"°- 5 m°- 4 , (24) 

R 

where we have suppressed the very weak dependence on 
y , to, and x. An upper, MRI-active layer can exist only if 
Z2/R < 1. If this layer does not exist, i.e., if to > 20a 1 - 2 
at our extremely simplified level of approximation, then 
the flow is in the soft state. 

4.3. Intermediate states 

If Z 2 < R and both lower and upper MRI-active zones 
coexist, then the flow is in an intermediate state, i.e., the 
spectrum has both soft and hard components (Fig. 3). 
Using eq. m to estimate the ratio of accretion rates, 
we have 

^ ~ 0.7y°- 3 a°- 9 m/ 0 - 7 , (25) 

TOl 

where we have suppressed the very weak dependence on 
m and x, but have retained the dependence on y for rea¬ 
sons that we will now explain. If m2 to 1, i.e., most 
power emerging from the cool equatorial zone, then hav¬ 
ing y ~ 1 in the upper layer would overcool it. Rather, 
we would expect y ~ m 2 / toi, implying 

rh 2 - 0.6a 1 ' 3 (26) 

is roughly independent of m 1 in this soft-intermediate 
state. According to our simple scalings, the soft- 
intermediate state exists for 0.6a 13 < to < 20a 12 , i.e., 
potentially over a large range of luminosities. 

When to 0.6a 1 ' 3 , to 2 is comparable to or larger than 
mi, y ~ 0(1) is required in order to cool the upper layer, 
and we have 

to 2 ~ 0.7a o - 9 TO5’ 3 . (27) 

Since the hard flux dominates over the soft flux, we iden¬ 
tify this state as the hard-intermediate state. 

We add the caveat that we have assumed r > 1 and 
subrelativistic electrons at z ~ R to obtain these ex¬ 
pressions. These assumptions may be violated, e.g., the 
optical depth in the upper radiating region may be 1, 


implying that y oc r rather than r 2 . Whether the elec¬ 
trons become mildly relativistic in this limit depends on 
whether y <C 1 (as in the soft-intermediate state) or 
~ 1. We note, however, that the hard spectral com¬ 
ponents in some XRBs in an intermediate state show 
much higher high-en ergy c utoffs than “pure” hard-state 
spectra (iMotta et al.ll2009h . 

4.4. Hard state 

When to < 0.6a 13 , to 2 rh\ and there is plenty of 
energy available to evaporate the cool disk and incorpo¬ 
rate its mass into the hot accretion flow. We therefore 
suppose that the cool disk vanishes, and seek a solution 
for a single-zone, magnetically supported flow in which 
radiation pressure is negligible com pared t o gas pressure. 
Equations (10)-(13) and (16) of iBeeelman fe Pringlel 
(I2007f ) remain valid, while p g is given by pc 2 . Assuming 
Coulomb coupling between electrons and ions as before, 
and cooling of electrons by unsaturated Comptonization 
with r > 1 and y ~ 0(1), we find that T\ T e for 
to > 10“ 2 p 1 / 2 ax, with 

^ ~ 0.5y- 3 / 7 a- 10 / 7 TO 8 / 7 x- 13 / 14 , (28) 

c 

(9 e ~ 0.04y 1/7 a“ 6/7 TO 2/ V /7 . (29) 

One can verify a posteriori that r > 1 for typical param¬ 
eters, justifying use of the subrelativistic version of the 
y— parameter. At lower values of to the electrons and 
ions would become thermally coupled, and thermaliza- 
tion of the radiation might occur. 

Once the flow is in the hard state, it need not return to 
an intermediate state when to > 0.6a 1 - 6 . As m increases, 
both Cg/c and 9 e are expected to increase. Dynamical 
constraints make it likely that c g /c will saturate at some 
value ~ 0(1), in which case the flow is essentially indis¬ 
tinguishable from a gas pressure supported flow. If 9 e 
also saturates at some value ~ 0(1), e.g., due to a steep 
increase in the cooling rate as the electrons become rel¬ 
ativistic, then we recover the usual criterion for main¬ 
taining a h ot (two-tem perature) accretion flow, m/a 2 < 
const. dRees et al.lll982f) . 

5. DISCUSSION AND CONCLUSIONS 

We have proposed that the spectral states in X-ray bi¬ 
naries arise as a consequence of the vertical structure of 
accretion disks whose dominant source of pressure sup¬ 
port is a toroidal magnetic field. As the field rises buoy¬ 
antly, its tension reaches a level comp ared to local gas 
press ure at which MRI is quenched (Pessa hfc Psaltisl 
1200511 . However, the Poynting flux from the equatorial 
dynamo continues to heat the overlying gas, the density 
of which declines so steeply that it undergoes a steep 
temperature inversion. This increase in gas sound speed 
can reignite MRI above the dead zone. The existence 
of one or two MRI-active layers, carrying independent 
mass fluxes and producing various ratios of soft and hard 
radiation, accounts naturally for the rich phenomenol¬ 
ogy of soft, hard and intermediate states. The strong 
toroidal magnetic field protects both the soft and inter- 
mediate states from e ither thermal or viscous instability 
(iBegelman fe Pringlell2007l ). 

Our model readily explains several features of XRB 
phenomenology that have proven hard to understand. 


















First, it explains the inevitability of soft-to-liard state 
transitions as the accretion rate declines. In models 
where both cold and hot accretion states can exist a t low 
accretion rates (|Chen et, al.1 119951 lEsin et, all 119971! , the 
conversion of the former into the latter virtually requires 
explosive heating to drive the density down by several 
orders of magnitude, so that electrons and ions can de¬ 
couple thermally. This transition is usually imposed by 
fiat. However, in our model the vertical stratification 
leads to a gradual transition, through soft- and hard- 
intermediate states, as the bulk of the energy generation 
passes smoothly from the cool layer to the hot layer with 
decreasing to. 

Once the cool zone disappears, the hard state consists 
of a single MRI-active layer, dominated by gas pressure 
and with thermal properties rather different than those of 
the hard-intermediate state. Because of the weaker den¬ 
sity stratification when the dead zone is absent, both the 
ion and electron temperatures are somewhat lower than 
in the intermediate states (although Tj can still greatly 
exceed T e ). Thus we can explain the curious observa¬ 
tional fact that the “pure” hard state, where no thermal 
component is apparent, often has a lower hig h-energy 
cutoff than either of the intermediate states (IMotta et al.1 
120091) . Furthermore, because the accretion flow in the 
hard state is fundamentally different from that in the 
intermediate state, its transition back to the intermedi¬ 
ate state during the rising phase of an outburst occurs 
at a different threshold. This can help to explain the 
hysteretic cycle observed in state transitions. 

The model discussed here predicts a local relationship 
between the mass accretion rate, a, and the state of the 
disk (specifically, how much energy is dissipated near the 
equator versus low density regions higher up). We re¬ 
quire that the disk pressure be dominated by toroidal 
magnetic field pressure at all heights, which in turn re¬ 
quires the presence of a poloidal field whose pressure 
exceeds about one percent of the gas plus radiation 
pressure. Observationally, the presence of jets during 
some phases of X-ray binary outbur sts str ongl y suggests 
that poloidal field is presen t (e.g. iBlandford fo Znaiekl 
119771 [McKinney et al.l [20l2h . but does not necessaril y 
imply that it be of large scale dParfrev et al.1 120151! . 
Disk winds, moreover, are inferred to be present in 
the soft and intermediate states of XRBs, consistent 
with a radially extended region that contains signifi¬ 
cant poloidal flux. Theoretically, maintaining even fairly 
weak poloidal fields in the inner disk is not necessarily 
straightforward; radial dragging of poloidal flux is in- 
efficient in geometrica lly thin disks (ILubow et al .1 H9941 : 
iGuilet fe Qgilviel [20121. and the extent to which an in¬ 
verse cascade can generate locally strong fields fro m zero 
net f l ux initial conditions is unknown (Sorathi a et, al.1 
120101: iBeckwith et al.1 [201 lh . We have previously sug¬ 
gested that the transition zone between geometrically 
thin and thick flows is a favora ble site f or sto chas- 
tic flux accumulation flBegelman fc Armitage!l20l4l! . and 
this mechanism may play a role in sustaining poloidal 
flux near the black hole. We note, in particular, that if 
the poloidal flux ever dropped below the level needed to 
sustain magnetic pressure dominance, the onset of ther¬ 
mal instability in radiation pressure dominated regions 
would create the geometry that is favorable for flux accu¬ 
mulation. More generally, changes in a driven by global 


changes in the poloidal flux appear necessary to explain 
why the intermediate state can exist at different levels of 
luminosity. 

A complete model for state transitions must also ac¬ 
count for the temporal properties of the states and their 
relationship to jet outflows. Strong rms fluctuations on 
short timescales, and lower frequency QPOs, are key 
features of the hard and intermediate states. In par¬ 
ticular, the hard and hard-intermediate states exhibit 
Type C QPOs, which are strong and relatively coher¬ 
ent (with Q —values as large as ~ 10), with a range 
of frequen cies (~ 0.1 — 10 Hz) correlated with sp ectral 
hardness dVignarca et al.l 120031 iStiele et ~ahl I2013T) and 
energy cutoff (E. Kalemci, private communication); while 
the soft-intermediate states sometimes exhibit less coher- 
ent QPOs with roughly fixed frequenc ies (^6 — 8 Hz) 
(ICasella et al.ll2005llBelloni et al.ll2011ll . The frequencies 
of these QPOs are all much lower than any dynamical or 
likely precession period in the region close to the black 
hole where most energy is generated, and we suggest that 
they reflect the quasi-periodic magnetic field reversals in 
the dynamo. Simulations indicate that the intervals be¬ 
tween flux reversals become quite long as the strength of 
the toroidal field increases, reaching tens of orbital peri¬ 
ods, if not longer, for magnetically dominated flows (Bai 
& Stone 2013; Salvesen et ah, in preparation). Such time 
scales are compatible with the observed QPOs, if they are 
produced at radii of ~ 10 — 100i? g . 

The high coherence of Type C QPOs suggests that 
they are either global modes, or are produced at a spe¬ 
cific point in the flow. We note, however, that for the ge¬ 
ometrically thick layers that likely modulate the QPOs, 
the inflow time scale is probably shorter than the field 
reversal time, implying that the polarity of the field is 
not produced locally, but rather propagates inward from 
some radius where the i nflow t i me an d field reversal time 
are comparable (lO’Neill et al.l [20Tlh . For example, this 
could be the radius where the flow switches from a thin 
disk to a thick, magnetically supported state. 

The model presented here is highly schematic. Our 
crude approximations include a very simple treatment of 
the vertical structure (taking a uniform magnetic pres¬ 
sure in each MRI-active zone, using piecewise power-law 
approximations for the pressures and densities, parame¬ 
terizing the critical question of how rapidly magnetic flux 
rises thorough the disk), and a minimal model for the 
electron-ion thermal coupling, which could differ from 
the Coulomb rate. The numerical values we have de¬ 
rived for the various thresholds are thus not expected to 
be accurate. We are encouraged, however, by the fact 
that given the reasonable assumption of a poloidal flux 
and resulting magnetic domination of the disk, the model 
shows promise for addressing a wide range of XRB phe¬ 
nomena. Moreover, the most uncertain aspects — the 
detailed vertical structure of strongly magnetized disks 
that include low density two-temperature regions, the 
accumulation and loss of poloidal flux in XRB disk ge¬ 
ometries, and the variability of emission in the hard state 
appear amenable to study using local and global sim¬ 
ulations that are feasible in the near-term. 
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APPENDIX 

EFFECT OF FIELD REVERSALS 


In many shearing-box simulations that show dynamo action, the direction of t he to ro idal he ld reverses qu asi- 
periodically, on a time scale which can be shorter than the rise time of the held ([Brandenburg et al .i f 19951 iStone et al.1 
Il99fit iDavis et~al]l20irt IQ ’Neill et al.ll20lit iSimon et alJl2012l : iBai fc Stonell2013ll.Tiang et al.ll2() 1 4aT . Regions of oppo¬ 
sitely directed toroidal held are separated by current sheets across which reconnection can take place, dissipating the 
magnetic energy more quickly than we have assumed in our analysis. This extra dissipation could lead to a steeper 
magnetic pressure gradient and heating concentrated more equatorially than in the unidirectional model. Here we 
present an extremely simple model that attempts to capture the magnitude and vertical distribution of this dissipation 
in a time-averaged sense. 

Suppose that the toroidal flux undergoes periodic reversals, with a duty cycle td ~ (£S2) _1 and £ ^ 1. These 
reversals, which originate near the equator, launch rising current sheets across which reconnection takes place with a 
speed 

Urec ~ Cl>A, (Al) 

where v\ is the Alfven speed on either side of the current sheet. The reconnection speed is defined as the speed with 
which the ambient magnetic field lines approach the current sheet from ei t her side, and recent s tudies of collisionless 
magnetic reconnection suggest that ( ~ 0.1 — 0.3 ( Zw eibel fe Ya mada 20 09tlUzdenskv et al.ll2010h . Thus, reconnection 
could occur rapidly compared to the rise time. As one current sheet sweeps by, the heat deposited per unit volume 
(and the magnetic energy density correspondingly lost) is 

AH ~ 2 ^ APb (A2) 

^rise 

where i; r i se is the rise speed of the current sheet, which need not be the same as the rise speed ub ~ yQ.z of the 
magnetic flux. Indeed, if the reconnection speed is required to be the same on both sides of the current sheet, then 
Uise cannot equal vb', moreover, there may be additional bulk flows of matter, both up and down (since the problem 
is now time-dependent) that should be included as the passage of the current sheet de posits momentum into the 

flow. (We suggest that these cyclic flows may be present in the density plot of Fig. 1 in lJiang et al.l[2014aT ). These 

complications are well beyond the simple level of our model, and to avoid the proliferation of undetermined parameters 
we will simply assume r) r j se ~ v rec , which then allows £ to cancel out of the problem. (We can subsume the uncertain 
validity of this assumption into the value of £.) 

The mean reconnection heating rate per unit volume, averaged over multiple cycles, is then 

H rec ~ 2,£Up B . (A3) 

Assuming that this energy comes from the Poynting flux generated by the MRI dynamo, we can replace eq. 0 by 

4-(pbVb) =ttvB[aBP+(as-2 Qpb\- (A4) 

az 

Using the same definitions for y and normalized quantities as before, we find that eq. © is unchanged, while the 
normalized magnetic energy flux equation 0 is replaced by 


2v~r(PBy) = («b - 2£)Pb + anPoP 
dy 


and the regularity condition becomes 

Defining a reconnection efficiency parameter, 

we find that eq. m becomes 


A-2a±fi-i. 

a b 

^2£ 

Q!B ’ 


vpb = 


1 + /3 0 - v 

The analytic approximation for pB , eq. m, is then replaced by 


(p-Pb). 


. l + Po-v ^-/VC+Ai-G 
Pb ~ ( 1 + 2(1 _|_ ^ o ) _ i,y 


(A5) 

(A6) 

(A7) 

(AS) 

(A9) 
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Because the Poynting flux now has an additional dissipative channel, the magnetic pressure decreases more steeply 
with height. The condition for magnetic pressure to dominate in the equatorial region (j3o < 1) is now p + £ < qb, 
but we find the same sensitive switch behavior as before, where a relatively small decrease in (77 + £,)/oib can trigger 
a drastic change in the structure of the disk-corona system, from gas+radiation pressure supported near the equator 
to overwhelmingly dominated by magnetic support everywhere. However, the condition for most of the energy flux to 
escape to z H is now more stringent, /3q < 1 — v. 

Among^ var ious s imulations extant which show a reversing accretion disk dynamo, we note that the models of 
iJiang et al.l (l2014al) show a much steeper drop in toroidal magnetic pressure, with distance from the disk core, than 
the pb oc z~ 2 estimated for our unidirectional field models with /3 0 1 (see their Fig. 2). We suggest that reconnection 

losses across the current sheets separating field reversals are responsible for this steepening. We note that their space- 
time diagram showing the “butterfly” pattern of reversing, buoyant fields (their Fig. 1, bottom panel) shows frequent 
reversals with one or more large-scale current sheets sweeping through the disk at all times. Their plots of volume 
dissipation rate as a function of column density, if fit by an alpha-model with T~L oc pb, allow us to estimate the 
reconnection efficiency parameter as v ~ 5 and 8 for models A and B, respectively. We note, however, that the 
reconnection in these calculations is driven by numerical resistivity, and does not necessarily reflect the degree of 
physical reconnection that would actually occur in such a disk. 

This approximate analysis is applicable to cases where field reversals occur on a shorter time scale than the rise time 
of the flux across a pressure scale height. Simulations suggest that this is the case for relatively weak toroidal fields 
that do not provide significant vertical suppo rt. For magne tically supported disks that are the main focus of this paper, 
field reversals appear to be very infrequent dBai fo Stond 120131 Salvesen et al., in preparation), and a unidirectional 
model as developed in Section 2 is appropriate. 
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